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Abstract 

A method to determine the full structure of the space of local operators of 
massive integrable field theories, based on the form factor bootstrap approach is 
presented. This method is applied to the integrable perturbations of the Ising 
conformal point. It is found that the content of local operators can be expressed 
in terms of fermionic sum representations of the characters xiQ) of the Virasoro 
irreducible representations of the minimal model M^^a,- The space of operators 
factorises into chiral components as Z = ^ but with the relation q = q^^. 
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1 Introduction 

In this letter we shall present a method in order to determine the operator content of an 
integrable massive field theory. It is based on the form factor bootstrap approach 0, 0], 
which provides a tool to classify all local operators of the theory under consideration. 
This is possible because the form factor equations are valid for every local operator of a 
theory. Therefore determining the whole solution space of these equations is equivalent 
to determine the content of local operators of the theory. 

As an example, we shall apply this method mainly to the integrable perturbations 
of the critical Ising model. The thermal perturbation can be described in terms of free 
massive Majorana fermions, or equally in terms of a bosonic particle which interacts 
through an S'-matrix S = —1. This theory consists of two sectors. One contains the 
monomials in the fermion fields, while the other sector contains fields which aquire a 
phase e*'^ when moved around the fermions. The basic fields in this second sector are the 
order and disorder fields a{x) and fi{x). 

The magnetic perturbation is described by a massive field theory containing 8 particles 
1^, and is intrinsically related to the algebra. There are no internal symmetries in 
this theory and all operators are expected to be mutually local. 

Let us summarise the description of an integrable theory in the form factor approach. 
We parametrise the momenta of the asymptotic states in terms of the mass m of the 
particles and the rapidity variables Pi 



Form factors are matrix elements of a local operator O between the vacuum and the set 
of asymptotic states. 



Their knowledge determines the correlation functions which can be expressed as 



= m cosh Pi , pI 



msinh/3j . 



J-f(A,/?2,...,/?n) = {0\O{0,0)\Z{P,),Z{P2),...,Z{P^)) 



(1) 



{Oiz,z)O{0)) 



(2) 




where r denotes the radial distance in the Euchdean space, i.e. r = ^Jx^ + xf, and s is 
the spin of the operator O. 

For systems with scalar particles, the functional equations known as the form-factor 
axioms are [jl], ||, ^: 



'' ei...eiei+i...e„l'^l' ■ ■ ■ i Hii • • • Hn) 

Se,e,+ Al3i - A+l) -^e^...e,+ ie,...e„(/5l, ■ ■ ■ , A, ■ • • Pn) 


(Fl) 


Sg„ denotes the spin of the particle e„; is the mutual locality index between 
O and Z^] w^-,,^ denotes the exchange properties of Z^. and Z^^. 




K..eM + A, . . . , /3„ + A) = e^^'J^Z.-eM^ • • - M 


(F3) 


T,, e (A + A, ...,/3i + A,/5i+i,...,/5„) =0(e^"'^l) for lAI ~ oo 

with max(5'^) < oo 


(F4) 


-I lim^,^^(/3' - + /3, A, . . . , = 
(l - e^-K+E-^.^ nr^, - A)) . . . , 


(F5) 


If the two-particle scattering amplitude exhibits a pole with the residue 

-^lim^,^,,..^X/?-<^.,)^e..,(/3) = {T%^ , then 

-I lim^,^^(/5' - /?)^,<^..,,,,^..„,„(/3i, . . . , /5' + zw^k, /? - • • • , /3n-i) = 

= r-,.^.^.......„(A,---,/3,---,/5n-i) 


(F6) 



Further we will use the following notation: We introduce the character x of 
conformal minimal model as 



where c and h denote the central charge and the conformal dimension of the primary op- 



erator respectively. Therefore x encodes the degeneracies in the Virasoro representation 
i, 

-| CO 

(Q')oo ^=_oo 

Further, define P(m, n) to be the number of partitions of n into numbers whose value do 
not exceed m. These partitions are generated by 

-1 CO 

-— = EP(m,n)g'^ , (5) 

where 

m 

iq)m ^ - 1') ■ (6) 

i=l 

2 The thermal perturbation of the Ising model 

The theory contains only one particle, and the form factors can be parametrised as ^, Q 

J^niPl ...(3n) = g„(Xi, . . . , X„)^n(a„)^^-»-^ f[ tauh ^A^M. , (7) 

where cij denote the elementary symmetric polynomials 0, Hn are normalisation con- 
stants and Xi = This parametrisation reduces the recursion relation (|l]) to 

Qn+2{-X,X,Xi,...,Xn) =0 , W = . 

The integer in the parametrisation introduces a grading into the space of 
operators. Increasing augments the divergence of the corresponding form factor, if the 
spin s of the operator O is kept fixed. Therefore, in analogy to CFT, one can define as 
the chiral operators those with the mildest ultraviolet behaviour, = We have 
therefore two discrete parameters in the massive theory, namely the integer and the 
spin s. The spin s can be related to the dimensions {h, h) of the conformal fields which 
we obtain taking the UV limit as usual through s = h — h. The parameter A^ has no 
obvious interpretation in the conformal limit. Its role will be clarified later on. 



2.1 Operators in the sector w = ^ 

The sector w = ^ contains the order and the disorder fields a and fi, and has been 
investigated in |§, |§] . The form-factors of both operators have been calculated and are 
given by = and Qn = const, in (0); further, in ^ the descendent operators of (t{x) 
have been determined. 

The method we want to describe does in general not resolve the recursion relations 
but rather counts the number of linearly independent solutions iteratively. Namely we 
use the fact that the number of solutions of the recursion relations at level n, i.e. for the 
form factor JF„, is given by the number of solutions at level n — 2 plus the dimension of 
the kernel of the recursion relations (|^). This amounts in a simple counting procedure 
which only involves the comparison of the degree of the polynomials Qn and that of the 
kernel. This has the advantage that one finds all solutions to the form factor equations, 
and not only those which can be considered as descendent operators of some non-trivial 
primary field. In the case of the Ising model the kernel of the recursion relation (H) is 
given by 

n 

JCn = YiiXi + Xj) . (9) 

i<j 

The total degree of this function is deg{lCn) = \n{n — 1). A Kernel solution is consistent 
for spin s only if deg{Qn) > deg{ICn), that is, Qn can be written as Qn = -R(crj)/C„, 
where R{ai) is an arbitrary combination of elementary symmetric polynomials of degree 
deg{Qn) - deg{K,n)- 

Using this fact and the generating function (P), one finds that the number of solutions 
at general spin s, is generated by 

„ir?i(m— 1) 

Fo^T. = & . (10) 

■m,odd \H>m 

where the exponents in the nominator corresponds to the spin value, where a new kernel 



solution enters. Note that this formalism leads to fermionic sum expressions [|10l of the 
character. 

Similarly one can analyse the even form-factors, being related to the disorder field ^. 
The counting procedure can be carried out analogously, and one finds that the number 



of chiral operators for spin s is generated by 

„im(m— 1) 

Go- E ^^T^-^if' ■ (11) 

m,even \H)m 

As it is well known, both order and disorder fields are related to the conformal operator 

As a next step we investigate solutions corresponding to higher N. Increasing we 
can gradually build up the whole space of operators in the massive theory. The number 
of operators at level A^ is generated by 

oo ^m{rn—l)—Nm 
m=l,odd 

in the odd sector, and by 



(9)r 



oo q^m(m~l)—Nm 

Gn= E 73 

m=0,even x'djm 

in the even sector. These functions satisfy the recursion relations 

Fn = Fn_i — (f Gn-\ , Gn = Gn~i + (f^ F^-i , (12) 

with the initial conditions Go = -^o = Xi,2{.<l)^ which can easily be solved to give F^ = 
Gn = nfcLi(l + (t)FQ. The whole space of states is generated by 

lim Fn = Xi,2(g)Xi,2(g) , (13) 

Af— >oo 

and similarly for Goo- 

Let us summarise the key features of the method described, which generalise also to 
other integrable massive field theories ||15[|: 



The content of chiral operators is given by fermionic sum expressions of the char- 
acters of the corresponding Virasoro irreducible representation (or of sums of the 
characters) . 

The whole space of operators is constructed by taking the limit N ^ oo. It can be 
formally decomposed into conformal characters as J2x{(l)x{(l) where the variables 
q and q are not independent as in CFT but q = q~^. 



Finally, one obtains for finite some 'finitized' expressions for the Virasoro char- 
acters. In this approach they appear naturally in the counting procedure of states. 
Note that the same expressions turn up in the study of corner transfer matrices 
where the finitization is related to the fact that one has a discrete system of finite 



size (see e.g. [|rT| for the thermal perturbation of the Ising model). It would be 



interesting to understand whether there is a deeper reason for this relation. 
2.2 Operators in the sector w = 

For the sector w = the kinematical recursion relation maps the form factors onto zero, 
^n{P + «7r,/3, . . .) = 0. It follows that form factors with different particle number n are 
not linked, and therefore any kernel solution will represent an acceptable form factor 
from the point of view of the form factor equations. 

We start our analysis with the space related to the even form factors. The fundamental 
operator is the energy density, whose form factor is given by JF2 = sinh |. Comparing this 
with the parametrisation (0), we find that in this sector the primary operator corresponds 
to a solution with = 1. 

Using the counting method as before, we find that the solutions at level N are gener- 
ated by 



n,even 



Expressing the lowest terms in Virasoro characters, one finds that /q = Xi,2{(l) and /i, 
which we conjectured to generate the chiral operators, is given by /i = Xi,i{q) + Xi,3{'l)y 
as expected. 

Since the form factors are not linked by the recursive equations for every operator in 
this sector there is only one contribution in the sum in (H), and the critical exponents can 
be calculated explicitly. One can show that all operators corresponding to solutions Qn in 
the space = 0, are operators which will scale in the ultraviolet limit as elements of the 
identity module, while the operators in the = 1 space form the space of descendents 
of the primary field 0i,3 f\. 

^The operators of the space iV = appear in the space iV = 1 as solutions Qn ~ cr„. 



In order to find the decomposition of tliis expression into cliaracters let us introduce 

(15) 



q 2 



n,odd 

Tliese generating functions satisfy the recursion identities 

9n = Qn-i + q'^ Jn-i , 



with initial conditions /o = Xi,i ^-^id Qq = Xi,3- Their solution is given by |T^, [H 
where 



fl= — OD 



-M(12At+l) 



fi=—oo 



:/i(l2^t+l) 



2N 
N -All 

2N 
N -1-4^ 



_ g:(3M+l)(4M+l) 



_ gr(3/.+l)(4M+2) 



2N 
N -l-Afi 

2N 
N-2-Afx 



(16) 



(17) 



:i8) 



and [ ] = denote the Gaussian polynomials. These expressions satisfy 

b 



(19) 



which means that for N ^ oo the space of local operators in this sector can be written 
as 

foo = Xi,i{q)Xi,i{<i) + Xi,3(g)Xi,3(g) • (20) 

Finally we analyse the odd form factors in the sector w = 0. From the parametrisation 
(^ we find that form-factors of the operators with the mildest ultraviolet behaviour 
(A^ = 0) will have half integer spin. We choose as our basic operator the fermion, whose 
form factor is given by J-'i = ai = . Carrying out a similar analysis as for the even 
form factors one finds that the generating function for the chiral operators is given by 



5^0 = Xi,3) where Gat is given in (|15|) . One can easily check that the ultraviolet dimensions 
are those of the corresponding conformal operators. 



The whole space of operators in this sector is generated by 

^oo = xi,i(^)xi,3(g) + gxi,3(g)xi,i(g) • (21) 

The space of operators has the same structure as the space of descendent operators of 
the conformal operators (|,0) and (0, |). 

Summarising, we found the descendent spaces of the fermions, the energy density 
and the identity operator in this sector. All are expressed in terms of the characters 
of the critical Ising model. It is interesting to note that the counting method selects 
automatically a base in the space of operators which is isomorphic to the corresponding 
Virasoro irreducible representation. 



3 The magnetic perturbation 

The perturbation of the critical Ising model by the operator 0i 2 couples the model to an 
external magnetic field at critical temperature. This system develops a finite correlation 
length and hence is massive. The on-shell structure of the theory has been determined 
in 0, and contains 8 scalar massive particles. The ^-matrix of the fundamental particle 
is given by 

SniPi - P2) = /2/3(/3i - /32)/2/5(/3i - P2)fiMPi - P2) 

with 

tanh ^(p — lan) 

All other S-matrix elements can be calculated by using the bootstrap equations. Because 
of the lack of space we are not able to give a detailed account of our investigations. We 
will present the major results, in order to show how the counting method can be applied 
also to relatively complicated systems, and give a complete presentation elsewhere [|15] . 



We start our analysis by applying the counting method to the fundamental particle 
and include the others by using the bound state equation (F6). In order to parametrise 
the form-factor conveniently, we introduce 

C(^, a) = exp { r dx coshx{l-a)sm^^^^^ ^ 
' /ox cosh I sinh x ' 



and 



< a >= sinh ^^f^ ~^ iira) sinh ~ iira) 



The form-factors corresponding to Ui particles of type 1 can then be parameterised 

as 

^„,(A,...,/?„J= (24) 

u n t 1 -^4- C(Aj; i)C(Ai) |)C(Ai; ]^) 

-'-'niVnA-^l, • • • , J-nJ 11 / n 3 i o 2 1 

<i i<j {XiXj)2 <1>2<|><|><^> 

where = f3i — j3j. Note that this parametrisation is an analytic continuation of the 



parametrisation of the form-factors of the particle of the BuUogh-Dodd model ||T6[. The 
degree of the kernel of the recursion relations is given by degiJCn^) = ^ni{ni — 1), while 
the degree of Qm is deg^Qm) = — 1), which leads to the g-sum expression 

2mi(mi-l) 

In order to determine the contributions to this g-sum expression arising from particle 2, 
we apply the bound state equation to (p4|). For the form factors with ni particles 1, and 
n2 particles 2 we find 

1 



^ 1 1,2. ..,2 — H, 



(n (r. r ))^^ (rr (r . r ))^^ 

^ Qni,n2 ("^l 5 • • • ) "^ni i -^ni+l ; • • • ; •^n\+n2 ) ^ (^6) 
^ Cil^iji ^^Cil^iji 5)C(Aj; 15) ^ 



,<j yXiXj)-2 < 1 >^< I >< I >< X > 



4 



MM(-.-.)^<t><i><^><i><i^> 

^ff^ <o>c(i)c(^)c(i)c(im^)c(^) 

f}, {^^^3f < 1 X I X I X ^ X 1^ X ^ X I X 1 X i X 1 > 
Comparing the dimension of the functions Q12 and Q22 with the dimension of the re- 
spective kernels, one finds that the exponent due to interaction of particles 1 and 2 is 
2m\ + Am\ + 4mim2 — 2mi — 4m2, modifying the expression (^) to 

/ 

E(mi,m2) I I I — 2r?ii— 4m2 
T^^^ ■ (27) 
{q)mMn.2 



4t 



Figure 1: The Eg Dynkin diagram, with the labelhng of the nodes used in the definition 
of the generating function 



Notice that the pole structure in ( PB] ) is not minimal (in the sense that they exhibit 
poles also in non-physical locations). This is no obstacle for the counting method since 
only the difference between the degrees of the functions Q and the kernels /C enters. 

Similarly it is possible to carry out the bootstrap for all particles. The corresponding 
expression of the form-factors becomes quite lengthly, and will be omitted. It will contain 
factors i = 1,...,8, corresponding to each of the 8 particles. The generating 

functions for levels Ni, . . . , can be found to be 

^MC-'^M-NC-^M 

Tni,N2,...,Ns = Tn = ^ -r~\ f~\ (28) 

where M denotes the multi- index M = mi, . . . mg and C^^ is the inverse Cartan matrix of 
the algebra Eg. The relevant labelling of the nodes of the Eg Dynkin diagram is depicted 
in fig. ^ The linear terms in the q sum expressions have been chosen in order to be 
consistent with the bootstrap. 

The lowest of these functions is % = Xi,i(Q') UHl [13]' corresponding to the most 
relevant operators, while the chiral operators are given by 



^0,1,0,0,0,0,0,0 = Xi,i{l) + XiAl) + Xi,3(g) • 
Similar higher functions can be expressed in terms of characters as for example 

^,0,1,0,0,0,0,0 = XiAq) + (1 + Q)Xi,2il) + XiAq) 

^0,0,0,1,0,0,0,0 = Xi,i(g) + (i + ^)xi,2(g) + (i + g)xi,3(g) 

To,o,o,o,i,o,o,o = (i + g')xi,i(g) + (i + g + g')xi,2(g) + (i + g + g')xi,3(g) • (29) 



These expressions correspond to the scalar partition function of the Ising model in its 
lowest orders in q. 

The generating functions (pH]) satisfy the recursion identities 

which allow to reconstruct the space of operators for higher levels. We have carried out 
extensive numerical calculations which show that for N ^ oo the functions 7/v approach 
the scalar partition function 

Z = Xi,i{q)Xi,i{Q) + XiM)XiM) + XiM)XiM) (30) 

as expected. It would be interesting to find exact solutions for these recursion relations, 
and determine their relation to the finitized characters for this model, which arise from 



carrying out the Bethe ansatz to the dilute As-model |17 . 



4 Conclusions 

We presented a method in order to determine the full content of local operators of an 
integrable massive field theory. We have discussed its application to the integrable per- 
turbations of the critical Ising model. It can be applied to any massive integrable theory. 



and in [|T5| we will discuss more examples explicitly. 

The method consists of a simple counting procedure and leads to fermionic sum 
expressions of the characters, or of sums of the characters. Further, by increasing the 
possible divergence of the form-factors, the whole space of operators can be obtained. 
We found that this space can be formally factorised into characters of the corresponding 
Virasoro irreducible representation of the form This is a quite remarkable 

result, since the chiral sectors are not independent in the massive model. This fact is 
reflected by the relation q = q~^. 

Though the structures of the conformal and massive theories are quite similar, there is 
one important difference. The grading in the space of conformal operators is introduced 
through the conformal weights of the operators, h and h, while in the massive model 
it is given through the spin s and the parameters A^. While the spin can be linked to 



the conformal weights as s = h — h the parameters lack such a direct interpretation 
in terms of the CFT weights. It is our behef that understanding this connection of 
to the Virasoro structure in the uhraviolet hmit, will give an indication of the algebraic 
structure determining the space of local operators in the massive theory. 
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